Instability of human societies as a result of conformity 
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We introduce a new model that mimics the strong and sudden effects induced by conformity in 
tightly interacting human societies. Such effects range from mere crowd phenomena to dramatic 
political turmoil. The model is a modified version of the fsing Hamiltonian. We have studied the 
properties of this Hamiltonian using both a Metropolis simulation and analytical derivations. Our 
study shows that increasing the value of the conformity parameter, results in a first order phase 
transition. As a result a majority of people begin honestly to support the idea that may contradict 
the moral principles of a normal human beings though each individual would support the moral 
principle without tight interaction with the society. Thus, above some critical level of conformity 
our society occurs to be instable with respect to ideas that might be doubtful. Our model includes, 
in a simplified way, human diversity with respect to loyalty to the moral principles. 
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INTRODUCTION 

It is commonly observed that strongly interacting 
groups of people may show very abrupt and unexpected 
behavior or opinion changes. These "transitions" may 
even lead people to behave against their inmost sense 
of good and evil. This is seen in many crowd phenom- 
ena, less harmfully, in fashion phenomena, or, more dra- 
matically, in the setting off of civil wars often leading 
to despotic regimes. Of course, these phenomena may 
have very different reasons, but we propose a very simple 
model that shows the same features and which may give 
some insights on the mechanisms of these regime rever- 
sals. Our hope is based upon similarities in the mani- 
festations of different non-democratic regimes which sug- 
gest that similar universal mechanisms may have played 
a role. 

All those regimes had powerful and merciless mecha- 
nisms of repression against people who were either ene- 
mies of the regime or merely neutral to it. However, one 
could argue that (see, for example, Q) the stability of 
those regimes was not based only on repression. In some 
cases, the honest belief of the majority of the population 
in the official ideas, or, at least the fact that people be- 
have as if it was the case, could play a major role. For 
example, as soon as this faithfulness became weak, the 
Soviet Union fell and the security forces were unable to 
stop this process. In many cases these ideologies drasti- 
cally contradict the moral principles of a normal human 
being. It is reasonable to think that those principles are 
shared by a vast majority of human beings. The experi- 
ence of the 20th century thus shows that under certain 
conditions these principles fail and the behavior of peo- 
ple seems to be controlled by some collective phenomena. 
Once more, we do not pretend that the model we have 



developed explains these historical events. Our modest 
goal is to show its analogies with human situations and 
to show how conformity can lead to unexpected extreme 
and paradoxical situations. 

These sociological or historical phenomena are often 
sudden, turmoil is drastic, and its apparent cause seems 
insignificant. For a physicist, this problem resembles a 
phase transition in a system consisting of interacting el- 
ements. In the system under study these elements are 
people and the most important part of their interac- 
tion is conformity. Conformity implies the ability, or the 
propensity, of people to change their opinion under influ- 
ence of the opinions of surrounding people. The confor- 
mity may be measured quantitatively by a simple exper- 
iment of the following type. A few pieces of paper each 
with its own number are demonstrated to a group of peo- 
ple. Each piece is either white or black. After the pieces 
are removed the experimentalist asks what is the color 
of piece #2. In fact, only one person in a group is being 
tested in this experiment; the rest of the group are also 
experimentalists, though the tested person is unaware of 
that. The tested person is asked the question after all the 
experimentalist's team presents the same wrong answer 
one after another. This experiment shows how difficult 
is to say that white is white after all other members of a 
group claimed that the white piece is black. The major- 
ity of people are unable to remain true to their conviction 
and will give a wrong answer (when the tested persons 
are alone, they almost always give the right answer) . As 
far as we know, these studies began in the middle 50th. 
Nowadays they flourish mostly in the purpose of mar- 
keting. Now everybody can take a free test to estimate 
roughly his/her level of conformity using Asch's Confor- 
mity Study Q. The role of conformity in a society has 
been considered previously 
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We do not wish to imply that conformity itself is some- 
thing bad. This is one of human properties which allows 
social life. However, we study below how conformity may 
lead to instabilities in human societies. 



STATISTICAL APPROACH 

Statistical physics is widely used for studying opin- 
ion dynamics in society in the framework of different 

models[! H0. 

There are good reviews of the field 8, 91. 
Our paper is devoted to a different aspect of the prob- 
lem, and we use here rather the Gibbs ensemble than 
dynamical equations. 

Our rough model considers the opinion of the society 
with respect to only one statement. We assume that all 
previous experience of the society including the inherited 
ideas of good and bad dictates a negative opinion with 
respect to that statement. However there is a confor- 
mity interaction between overlapping groups of the soci- 
ety and also there are some individuals who, in spite of 
everything, have a positive opinion with respect to the 
statement. Their positive opinion is a little bit stronger 
than the negative opinion of other people. Following the 
terminology of Mobilia H we shall call these individu- 
als "zealots". The number of zealots is not fixed and is 
determined from the thermal equilibrium. We interpret 
the phase transition in this system as social instability. 
At small conformity the majority has a negative opin- 
ion while at high conformity the opinion of the majority 
becomes strongly positive, i.e. the majority of people 
become zealots. 

The Hamiltonian of the model reads: 

H = ^2 4>i qi - - ^2 V hi lilj ■ W 

i i^j 

For simplicity it is formulated on a two-dimensional 
lattice. Variables qi may take values —1/2 or 1/2 + a. 
They represent the opinion of the person on a site i with 
respect to a given statement. In terms of the previous 
terminology zealots have positive opinions and value a > 
0. This parameter describes some extra strength of the 
positive opinion, the negative opinion being the " normal" 
one (whether because it corresponds to the true color 
of the piece of paper, or to the opinion dictated by the 
"universal" notion of good and evil. . . ). In the absence of 
conformity (V = 0), at low temperatures all the opinions 
are negative (within thermal fluctuations). 

The random values <pi are positive with a Gaussian 
distribution F(<f>) = V2/(^A)e^q>-(<j) 2 /2A 2 ). They 
describe the diversity of the individuals with respect to 
the statement under consideration. Indeed, considering 
only the first term in Eq. (T|. one finds that the energy 
H has a minimum when all values of qi are negative. The 



second term describes interaction and conformity. It has 
a minimum when all qi have the same sign, in other words 
when all people think in the same way. Moreover, at any 
positive value of a, it prefers all people to be zealots. For 
simplicity we consider only the interactions with the four 
nearest neighbors on a square lattice and Vij = V. 

Eq. (JIJ can be represented as an Hamiltonian of the 
Ising model in some fictitious magnetic field[lfj| . Intro- 
ducing the new variables Si = qi — a/2 one gets: 

ff =E*-^ Wi+c ' (2) 

i ijij 

where: 

j 

is a random "magnetic field" at a site i, C is is irrele- 
vant constant, and the values of Si are ±1/2(1 + a). 

FIRST ORDER PHASE TRANSITION 

We now show that our system has a first order phase 
transition. At zero temperature it can be derived exactly. 
To find the transition point we compare the energies of 
the two states. In the state I all qi = —1/2. It has the 
energy per site: 



& = -jW-y (4) 

In the state II all = 1/2 + a which means all people 
are zealots. It has an energy per site: 

^=(!+a)<*>-£(i+a) a . (5) 

Here (</>) = /" 4> F{4>) ckf> and z = 4 is the number 
of nearest neighbors taking part in the interaction with 
a given site. When calculating the interaction per site 
we must take V/2 since this contribution belongs to two 
sites. It is easily found that: 




One can see that, at small interaction az K <C 1, the 
state I has the lower energy while in the opposite case the 
state II has the lower energy and represents the ground 
state of the system. This is a typical situation for a first 
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order phase transition. It happens at Ej = En. Using 
Eqs. 1415161) one finds that the transition occurs at: 



1.6 

a z 



(7) 



and the energy in the transition point is given by the 
equation: 
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where N is the number of sites in the system. Note that 
this result can be obtained from the condition < h >= 0, 
where h is given by Eq.J3J). 

If conformity is small (-j < (^) t ), the total energy 
can be found from Eqs. 1 ll(il) . It has the form: 
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The total energy at large conformity can be found from 
Eqs. One obtains: 



H H = VN 
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MONTE-CARLO SIMULATION AT FINITE 
TEMPERATURE 

A Metropolis code has been written for the simu- 
lation of the Hamiltonian (Eq. on a 2D (L x L = N) 
square periodical lattice. At each Monte-Carlo step, a 
site k of the lattice is chosen at random. The value 
of qk of the site is flipped with a probability Pk — 
Min(l, exp —SHk/k^T) where SHk = tk 5<lk is the energy 
cost of flipping. Here: 



£fc = 4>k 
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(11) 



It is well known that, after an initial number of Monte- 
Carlo steps that are necessary for relaxation, the averag- 
ing of any function of the qi over the subsequent set of 
distribution of qi, obtained by this way, is equivalent to 
averaging over the Gibbs ensemble. Depending on the 
conditions, the operation is performed from one hundred 
to one million times per site. The simulation has been 
run for different values of the independent parameters 
T, -j, a, starting from different initial distributions of 
the qi. In what follows we use, instead of T, a dimen- 
sionless temperature, so that T — > ksT/V . 



The main results of the simulation are shown in Fig. 
^ The time averaged value (q) is plotted against ^ at 
various values of the temperature at a = 0.1. One can see 
that, at low temperatures and low conformity the point 
of view of the population with respect to the statement 
under study is negative. At strong conformity {Xr — > oo), 
all members of the population become zealots. At low 
temperatures one can see a wide hysteresis. At higher 
temperatures it disappears and all curves intersect at one 
point at -j ~ 4. This value is in a good agreement with 
Eq. at z — 4 and a = 0.1. The value of (q) in the 
transition point is w 0.05 (i.e. (-1/2 + 1/2 + a)/2). The 
hysteresis loop collapses with increasing temperature. 




FIG. 1: Average public opinion as a function of ^ at differ- 
ent temperatures (circles: 0.4, squares: 0.5, triangles: 0.6, 
plus: 0.7, crosses: 0.8, stars: 0.9, diamonds: 1.0) as given by 
the simulation. The results obtained for increasing confor- 
mity are connected by a solid line, those corresponding to the 
backward path are connected by dotted lines. For tempera- 
tures greater than 0.8 both paths are superimposed (within 
numerical fluctuations) . 



Fig. |3shows the simulation results for (^) t as a func- 
tion of a at different temperatures. One can see that, at 
the transition, it reproduces very well Eq. J7J which was 
obtained at zero temperature (the same result can be ob- 
tained writing that at the transition (e) =0, Eq. 111(1 ). 
The simulation shows that Eq. (0 is applicable over a 
wide temperature range. This follows also from the fact 
that curves corresponding to different temperature cross 
at one point in Fig. ^ 

However the energy at the transition point differs sub- 
stantially from the zero-temperature energy given by Eq. 
(JHJ. The transition energy H t is shown in Fig. To 
avoid the auto-magnetisation effect at low temperature, 
the initial directions of the spins are sorted at random. 

It is also interesting to compare the behavior of the 
energy as a function of K at a finite temperature with 
the zero-temperature behavior as given by Eqs. I|9I10|I . 
These results are given in Fig. 0] 
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FIG. 2: Transition value (-j) t as a function of a. The line 
corresponds to Eq. J7J and the dots to the simulation results 
obtained for different values of the temperature (T € [0, 1]). 
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FIG. 4: Total energy as a function of ¥ in the case T = 
0.7, a = 0.1 obtained on a 100 x 100 lattice (open dots). The 
filled dots correspond to one quarter of the energy obtained 
for a 200 x 200 lattice. The zero-temperature results given by 
Eqs. H9I101 are shown by a solid line. 
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FIG. 3: Transition value Ht as a function of a. The line corre- 
sponds to Eq. derived at T = 0. The points are obtained 
from the simulation for different values of the temperature 
(diamonds: T = 1, squares: T = 0.7, circles: T = 0.1). 



DISCUSSION OF COMPUTATIONAL RESULTS 

As far as we know, the Hamiltonian given by Eq. 
has never been studied before. At A = a = it coin- 
cides with the ferromagnetic Ising Hamiltonian and it has 
a second-order phase transition in the two-dimensional 
case. We claim that our Hamiltonian has a first order 
transition at a point given by Eq. J7J). It is obvious at 
T = 0. Convincing argument is a singularity in the be- 
havior of energy as a function of Xr shown in Fig. 0]that 
we interpret as a discontinuity of the first derivative of 
the energy in the transition point. Our understanding 
of the hysteresis is that our modeling in real time goes 
very slowly and therefore we may see a continuation of 



the phase from which we start above the transition point. 
That is typical for the first order transition and reminds 
super-cooled liquid. 

On the other hand, we study a very low value of a 
and transition occurs at small value of A. Therefore one 
should expect that it should be close to the second order 
phase transition in the Ising model. Therefore we may 
observe clusters of zealots near the phase transition. As 
an example the distribution of the zealots at the thresh- 
old of the transition is shown in Fig. [5J One can see that 
the system has a large correlation radius which is typical 
for the second order phase transition. Thus, we think 
that we have a weak first order phase transition. 



CONCLUSIONS 

We propose a simple Hamiltonian that models the 
drastic opinion changes that can be experienced by hu- 
man groups submitted to strong mutual influence, even 
changes that contradict a cultural achievements of the 
past. This model shows how, due to the conformity, the 
group may express an opinion opposite to the opinion 
which each of its isolated members would have. The 
model contains the following parameters: T, A, V, a, z 
that can be determined by sociological methods. The 
level of conformity may be checked by experiments sim- 
ilar to those described in the Introduction. The "tem- 
perature" of the society could be found by studying time 
fluctuations of the public opinion. 

Our model is oversimplified, and the Eq. (JJJ) for the 
critical value at the first order transition point may be 
not accurate. Nevertheless, if the mechanism of collec- 
tive phenomenon in the human society can basically be 
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FIG. 5: Distribution of the zealots at T = 0.62, \ = 
, N = 10 4 , a = 0.1. The final distribution of the 
zealots are indicated by black squares (the 100x100 lattice 
is periodical). 



described our model, the ratio (-j) t as given by Eq. Q 
might be an important characteristic of group phenom- 



ena. 
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